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ABSTRACT 

We have measured the positions of 74 pulsars from regular timing observations 
using the Nanshan radio telescope at Urumqi Observatory between 2000 January and 
2004 August (MJD 51500 - 53240). Proper motions were determined for these pulsars 
by comparing their current positions with positions given in pulsar catalogues. We 
compare our results to earlier measurements in the literature and show that, in gen- 
eral, the values agree. New or improved proper motions are obtained for 16 pulsars. 
The effect of period fluctuations and other timing noise on the determination of pul- 
sar positions is investigated. For our sample, the mean and rms transverse velocities 
are 443 and 224 km s" 1 respectively, agreeing with previous work even though we 
determine distances using the new NE2001 electron density model. 

Key words: astrometry - pulsars: general - methods: data analysis 



1 INTRODUCTION 



Most pulsar proper motions have been obtained using in 
terferometric observations l|Lyne, Anderson & Salter 1982 



Harrison, Lyne fc Anderson 1993| lFbmatonT^eT^L*T997 
Bris ken^tHu^OOSlT However, pulsar timing also al- 



lows measurement of very accurate pulsar positions and 
hence proper motions ^Ma nchester , Taylor fc Van 1974] 
IToscano et al. 1 999). This method was considered to 
produce inaccurate results for young pulsars whose timing 
residuals may be dominated by glitches and timing noise. 
Recently, Hobbs et al. (2004) showed that with data 
spanning more than six years it is possible to model the 
timing noise by fitting harmonically related sinusoids. The 
spectrum of this timing noise has little power at periods 
corresponding to one year and so any timing noise can be 
removed from the timing residuals whilst leaving the signa- 
tures of an incorrect position determination or unmodelled 
proper motion. In this paper, we use the simple method of 
determining pulsar proper motions by comparing current 
positions obtained with nearly five years of timing data 
from the Nanshan radio telescope at Urumqi Observatory 
in China to earlier catalogued values. Having only five years 
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of observations means that it is possible to fit just a few 
harmonically related sinusoids to the data, which in many 
cases leaves significant structures in the timing residuals 
after the whitening procedure. We therefore choose to 
whiten the residuals using polynomials instead of sinusoids. 
We compare the results from this analysis with published 
measurements of the proper motions and consider the 
effects of timing noise in these relatively short data spans. 



Proper motion measurements have shown that pulsars 
are high-velocity celestial objects with typical velocities of 
several hundred km s _1 j-Lyne fc Lorime" 1994), but it is 
not clear how the pulsars achieve such large velocities. The 
exact form of this velocity distribution is also not well un- 
derstood with Arzoumanian, Chernoff & Cordes (2002) find- 
ing evidence for a bimodal distribution, but Hobbs et al. 
(2005) finding no such evidence. An accurate determination 
of the proper motions is therefore essential in understand- 
ing the velocity distribution. However, the literature con- 
tains many examples of inconsistent proper motion and po- 
sition measurements. For instance Lyne, Anderson & Salter 
(1982) obtained fi a = — 102(5) mas yr -1 for the proper mo- 
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tion in right ascension 1 for PSR B1133+16 whereas Brisken 
(2001) published that fj, a = -74.0(4) masyr" 1 . Similarly 
Hobbs et al. (2004) determined a proper motion of fi a = 
— 10(12) masyr for B1552— 31 that is inconsistent with the 
value of 61(19) masyr -1 given by Brisken et al. (2002). It is 
therefore essential to provide many checks on the correct- 
ness of existing proper motion values as well as providing 
new or improved values. 

In §2, we briefly describe the observing system at 
the Nanshan radio telescope. In §3 we present our analysis 
method and the derived positions and proper motion val- 
ues. The discussion in §4 is divided into i) an analysis of 
the effect of timing noise on position determinations, ii) a 
comparison of our proper motion results with earlier pub- 
lished values, iii) the velocities obtained using our complete 
sample of proper motions and iv) new and improved proper 
motions from our work. 



2 OBSERVATIONS AND ANALYSIS 

Timing observations of 74 pulsars using the Nanshan 25-m 
radio telescope operated by the Urumqi Observatory, Na- 
tional Astronomical Observatories of China, have been reg- 
ularly carried out from 2000 January to 2002 June using a 
room temperature receiver at 1540 MHz. Since 2002 July, 
284 pulsars have been monitored approximately once every 
nine days using a dual-channel cryogenic system that re- 
ceives orthogonal linear polarisations at the central observ- 
ing frequency of 1540 MHz. The pulsars were chosen to be 
observable from the Observatory (which has a lower declina- 
tion limit of —43°) and to have mean flux densities greater 
than 0.5 mjy. In this paper, we consider results for the 74 
pulsars with mean flux densities larger than 2 mjy and dis- 
tances within 8.5 kpc over the data span 2000 January to 
2004 August. 

After down-conversion, the two polarisations are each 
fed to a niterbank consisting of 128 contiguous channels 
each of width 2.5 MHz. The signals from each channel are 
square-law detected, filtered and one-bit sampled at 1 ms 
intervals using a data acquisition system based on a PC 
operating under Windows NT (full details of the observ- 
ing system are provided in Wang et al. 2001). The data 
are folded synchronously at the pulsar topocentric period 
for sub-integration times of between one and four minutes 
and written to disk. Time-tagging of the pulse integrations 
is provided by a hydrogen maser clock calibrated using the 
Global Positioning System (GPS) signals and a latched mi- 
crosecond counter. Off-line programs sum the orthogonal 
polarisations and dedisperse the data into eight sub-bands. 
Observations of a given pulsar typically consist of four sub- 
integrations with observation lengths of between 4 and 16 
minutes depending on the flux density of the pulsar. The 
folded profiles are stored on disk for subsequent processing 
and analysis. 

We use the psrchive software package 
(Hota ~van Straten fc Manchester 2004| | to obtain the 
pulse times of arrival (TOAs) from each observation. To 

1 Uncertainties given in parentheses after a quantity refer to the 
last quoted digit and are la values. 



do this, the pulse profiles from each integration are cross- 
correlated with high signal-to-noise templates to produce 
topocentric TOAs. These are subsequently processed using 
the standard timing program tempo 2 with the improved 
Jet Propulsion Laboratory ephemeris DE405 3 . TOAs are 
weighted by the inverse square of their uncertainty, with 
these being scaled (if necessary) so that the fitting routine 
gives reduced \ 2 values close to 1.0. Uncertainties in the 
fitted parameters for most pulsars are taken to be twice the 
standard errors obtained from tempo. However, for pulsars 
whose timing residuals are dominated by timing noise the 
uncertainties in the fitted positions are taken to be three 
times the formal tempo errors. These error estimates are 
discussed in more detail in §4. 



3 RESULTS 

In this paper we present the results from an analysis of TOA 
data for 74 pulsars. These results include new positions from 
which we determine proper motions and velocities. For each 
of the 74 pulsars, the position was measured by setting the 
reference epoch of the timing model to an integral MJD near 
the centre of the data span. We fit the pulsar position, ro- 
tational frequency and its first derivative. For pulsars whose 
timing residuals are dominated by timing noise we fit multi- 
ple frequency derivatives to model and hence effectively re- 
move the timing noise. For instance, FigureQplots the post- 
fit timing residuals for PSRs J0139+5814 and J1709-1640 
after fitting for the rotational frequency and its first deriva- 
tive and after fitting multiple frequency derivatives. Among 
the 74 pulsars, the residuals of 26 are dominated by timing 
noise and ten glitches in total were found in five pulsars. 
Some timing results were described by Zou et al. (2004) and 
others will be reported in a subsequent paper. The typical 
post-fit residuals of the 74 pulsars are several hundred [is. 
All of the pulsars in our sample are solitary and none are 
millisecond pulsars. Most have ages in the range 10 4 — 10 7 yr 
and all are within 9 kpc of the Sun. 

In order to obtain proper motions we compared our 
results with early previously published position estimates. 
However, the very earliest position estimates often used rela- 
tively poor instrumentation and so can have large uncertain- 
ties. For each pulsar, we therefore select the earliest avail- 
able position estimate that has an uncertainty less than 3". 
For the majority of our pulsars, we obtained these positions 
from the Taylor, Manchester & Lyne (1993) catalogue. For 
PSRs B0525+21, B1749-28, B1754-24 and B1821-19 we 
use positions given by Hobbs et al. (2004). For most pul- 
sars, the epoch of the positions given by the catalogue is 
more than ten years earlier than the epoch for our data. 
The maximum time span is 33 yr. 

Table 1 gives the pulsar names in J2000 coordinates 
and B1950 coordinates in the first two columns. The follow- 
ing five columns give, respectively, the epoch of the position 
measurement, the right ascension and declination from the 
timing analysis, the number of TOAs in the nearly five-year 
data set and the time span between the two position de- 
terminations. The final four columns refer to the previously 

2 See http://www.atnf.csiro.au/research/pulsar/tempo 

3 See http://ssd.jpl.nasa.gov/iau-comm4/relateds.html 
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PSR J0139+5B14 



PSH J1709-1640 




published position and give the epoch, right ascension, dec- 
lination and a reference for the measurement. 

Proper motions derived by comparison of the Nanshan 
timing results with previously published positions for the 74 
pulsars as well as the most reliable previously published tim- 
ing and interferometric proper motions for these pulsars are 
given in Table 2. For most of the pulsars in our sample we 
determine the pulsar distance D, given in the third column 
of the table, using the dispersion measure and the Cordes 
& Lazio (2002; hereafter NE2001) electron density model. 
The exceptions are 13 pulsars for which interferometric an- 
nual parallax measurements are available 4 . The next three 
columns give the proper motion in the right-ascension direc- 
tion (Ha) the declination direction and the total proper 
motion (/Xtot) derived from the data listed in Table 1. If 
proper motions already exist in the literature, we also pro- 
vide the most reliable proper motions available from timing 
(superscript T) and from interferometry (superscript I). All 
of the timing measurements are from Hobbs et al. (2004) 
except that for PSR J0738— 4042 which is from Siegman et 
al. (1993). References for the interferometric proper motions 
are given in the 11th column. Offsets of fi c a and fig from the 
most accurate of the previously published proper motions 
normalised by the uncertainty in fi% and fig, respectively, 
are given in the next two columns (e a and eg). In the final 
column, values of the pulsar transverse velocity derived from 
Utot and the distance. 



4 DISCUSSION 

4.1 Errors in position measurements based on 
pulsar timing 

The precision and reliability of proper motions determined 
by comparison of timing positions with other positions are 

4 ATNF Pulsar Catalogue 

|http: //www.atnf.csiro.au/research/pulsar/pgrcat| 



clearly affected by the presence of timing noise. In this sec- 
tion we attempt to quantify the uncertainty in positions de- 
rived from pulsar timing as a function of the amplitude and 
spectrum of the timing noise and the length of the data span. 
Timing noise consists of two components: a flat-spectrum or 
'white' component, usually dominated by receiver noise, and 
a steep-spectrum or 'red' component normally dominated 
by fluctuations in the pulsar spin rate. The red or intrinsic 
component is generally larger in young pulsars and is often 
quantified by the stability parameter as 

A 8 = log(i-|^|t 3 ), (1) 

where v — 1/P and v is its second time derivative and 
t — 10 8 s l Arzoumania nTt al. 199H . These authors show 
that Ag is approximately proportional to P 1 ^ 2 , where P 
is the first time-derivative of the pulsar period. Timing 
noise can also be characterised by the power spectrum of 
the timing residuals S(y) ~ v a , where a — —2 if the 
noise results from a random walk in pulse phase, a = 
—4 for frequency noise and a — — 6 for slowing down 
noise (Boynton et al. 1972; Groth 1975). However, more 
recent work ( |D'Alessandro, D cshpand e fc McCulloch 1997^ 
has shown that the power-law spectra are often complex 
with several power-law segments. 

We investigate the effect of timing noise on position 
determinations by simulating pulse TO As containing a mix- 
ture of white and red noise using the tempo2 software pack- 
age jHobbs, Edwards fc M anches ter 2005^ . The red noise 
has a spectral index a = — 3 chosen to give an approxi- 
mate representation of the timing noise seen in our residuals 
(a = for white noise). TOAs were computed for data spans 
(T) of 2, 3, 5, 10 and 20 years at 10-day intervals and ad- 
justed so that the residuals matched the desired noise model. 
This process was repeated to give 100 sets of TOAs for each 
of the five data spans. Figure shows the distribution of 
offsets of the right ascension obtained from the timing fit 
compared to the nominal value for data sets dominated by 
white noise (a) and red noise (b). Figure^ shows that for 
white noise the error in the timing position decreases with 
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Table 1. Positions of 74 pulsars from Nanshan timing observations and previously published positions. 
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* Quoted errors are three times the formal TEMPO values (see text); others are two times the formal TEMPO values. 

References for the previously published data arc: (1) Taylor, Manchester & Lyno, 1993; (2) Arzoumanian ct al., 1994; (3) Johnston ct al., 1995; (4) 
Manchester ct al., 1996. (5) Hobbs et al., 2004. 



increasing data span, approximately proportional to T~ ' 
as expected. Furthermore, the standard deviation given by 
tempo accurately represents the uncertainty, with 68% of 
offsets within ±1<t. Figure |5Jj shows that, for red noise, 
the tempo standard deviation greatly under-estimates the 
actual scatter, especially for short data spans, and that it 



doesn't decrease as expected with increasing data span. For 
a timespan of 5 years, only 27% of the offsets are within ±1<t 
of the TEMPO solution. 

For the five-year data span relevant to the present ob- 
servations, we plot in Figure |31 an error factor, defined as 
the rms scatter of the points in Figure divided by the 
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Table 2. Comparison of pulsar proper motions derived from Nanshan timing results with previously published values. Proper motions 
are given as mas yr -1 on the sky, i.e., fi a = a cos S. 



PSR. J 


PSR B 


Dist 






M tot 


fa 


T 


fa 


fs' 


Rcf . 






5 








(kpc) 


(mas yr ) 


(mas yr 1 ) 


(mas yr" 1 ) 


(mas yr - 1 ) 


(mas yr - 1 ) 


(mas yr" 1 ) 


(mas yr" 1 ) 










(km a" 1 ) 


0034-0721 


0031-07 





4 


7(27) 


18(56) 


19(53) 


-16(26) 


17(53) 


-102(74) 


-105(78) 


6 


0.2 





,0 


38(100) 


0139 + 5814 


0136+57 


2 


.9 


— 21 (4) 


— 5(5) 


21(4) 






— 11(5) 


— 19(5) 


1 


1. 1 


1 


4 


290(52) 


01414-6009 


0138+59 


2 


.2 


-29(34) 


-12(23) 


32(33) 


10(14) 


-5(15) 








0.8 





.2 


330(340) 


03324-5434 


0329+54 


1 


.1 


12(3) 


5(4) 


12(3) 


12(4) 


-14(7) 


17.0(3) 


-9.5(4) 


8 


1.5 


3 


.3 


63(15) 


0358 + 5413 


0355 + 54 


1 


.1 


1(5) 


19(7) 


19(7) 


8(3) 


19(6) 


9.2(2) 


8.2(4) 


9 


1.6 


1 


.4 


99(37) 


0452-1759 


0450-18 


2 


.4 


1(4) 


8(7) 


8(6) 


10(4) 


3(6) 


12(8) 


18(15) 


3 


1.1 





.4 


94(70) 


0454 + 5543 


0450+55 


o 


.7 




14(9) 






— 13(12) 




— i7(2) 










210(10) 


0528 + 2200 


0525 + 21 


1 


6 


-24(20) 


-520(570) 


520(570) 


-38(15) 


200(300) 


-20(19) 


-7(9) 


1 


0.1 





.9 


5200(2800) 


0534 + 2200 


0531 + 21 


1 


.7 


-75(44) 


570(590) 


580(580) 


— 


— 


-16(11) 


-2(8) 


4 


1.1 


1 


.0 


4800(4800) 


0543 + 2329 


0540+23 


2 


.1 


6(42) 


140(1100) 


140(1100) 


22(32) 


-400(1100) 


19(7) 


12(8) 


1 


0.3 





.1 


1400(11000) 


0612+3721 


0609+37 


o 


9 


14(5) 


10(17) 


17(11) 


5(5) 


6(22) 








0.9 





1 


71(44) 


0630 — 2834 


0528 28 


1 


.5 


— 46(5) 


23(7) 


51(5) 


— M (7\ 


17(8) 


44,6(9) 


20(2) 












0738-4042 


0736-40 


2 


6 


-48(2) 


35(2) 


59(2) 


-60(10) 


40(10) 


-14.0(12) 


13(2) 


7 


11.9 


5 


.5 


740(18) 


0742-2822 


0740-28 


2 


.1 


-14(9) 


23(12) 


27(11) 


18(20) 


-49(21) 


-29(2) 


4(2) 


3 


1.3 


1 


.3 


260(110) 


0814+7429 


0809+74 





.4 


31(13) 


-52(19) 


61(18) 


— 


— 


24.02(9) 


-44.0(4) 


8 


0.5 





.4 


120(36) 


0820-1350 


0818-13 


2 





18(9) 


-44(9) 


47(9) 


21(3) 


-41(5) 


9(9) 


-47(6) 


3 


0.3 





.2 


440(85) 


0826 + 2637 


0823 + 26 





.4 


67(4) 


— 82(14) 


110(11) 


65(3) 


— 108(8) 


61(3) 


— 90(2) 


4 


0.9 


o 


,5 


180(18) 


0837-4135 


0835-41 


1 


1 


7(5) 


-6(5) 


9(5) 






-2.3(18) 


-18(3) 


2 


1.4 


1 


.5 


44(24) 


0846-3533 


0844-35 





6 


93(72) 


-15(65) 


94(72) 


54(70) 


-8(85) 


— 


— 


— 


0.3 





.0 


250(190) 


0922+0638 


0919+06 


1 


2 


21(14) 


28(34) 


34(31) 


— 


— 


18.8(9) 


86.4(7) 


7 


0.1 


1 


.6 


200(170) 


0953+0755 


0950+08 





3 


-22(14) 


-39(41) 


45(36) 


-3(3) 


26(7) 


-2.09(8) 


29.46(7) 


8 


1.4 


1 


.7 


55(45) 


1 136 + 1551 


1133+16 





.4 


— 61(11) 


380( 10) 


380( 10) 


— 69(3) 


327(5) 


— 74.0(4) 


368. 1(3) 


8 


1.2 


1 


o 


650(17) 


1239 + 2453 


1237+25 





9 


-100(7) 


65(10) 


120(8) 


-108.8(14) 


50(3) 


-104.5(11) 


49.4(14) 


7 


0.4 


1 


.4 


490(34) 


1509 + 5531 


1508+55 


1 


.0 


-75(9) 


-49(9) 


89(9) 


-86(12) 


-65(11) 


-70.6(16) 


-68.8(12) 


7 


0.4 


1 


.9 


420(42) 


1645-0317 


1642-03 


1 


.1 


-4(3) 


10(8) 


11(7) 


— 


— 


-3.7(15) 


30.5(16) 


2 


0.1 


2 


.2 


57(29) 


1703-3241 


1700-32 


2 


.3 


170(130) 


-240(300) 


290(250) 


-26(17) 


62(104) 


— 


— 


— 


1.4 





.7 


3200(2800) 


1705— 1906 


1702 — 19 


o 


.9 


— 77(14) 


— 15(120) 


78(26) 


— 78(9) 


— 1 16(82) 








0.0 


o 




330(110) 


1707—4053 


1703-40 


4 


5 


3(21) 


35(61) 


35(61) 






— 


— 


— 








760(1300) 


1709-1640 


1706-16 





.8 


-9(4) 


13(34) 


16(28) 


— 


— 


3(9) 


0(14) 


3 


0.9 





.3 


64(77) 


1721-3532 


1718-35 


5 


6 


-43(24) 


-270(93) 


270(91) 


-6(105) 


-600(600) 


— 


— 


— 


0.3 





.5 


7200(2400) 


1722 — 3207 


1718-32 


2 


.4 


-1(5) 


-40(27) 


40(27) 


-5(8) 


5(54) 


— 


— 


— 


0.3 





.6 


450(300) 


1740 — 3015 


1737 — 30 


2 


.7 


— 170( 190) 


55(380) 


180(210) 


















240 0(2700) 


1741-3927 


1737-39 


3 


.2 


20(15) 


-6(59) 


21(22) 


— 


— 


— 


— 


— 








320(340) 


1745-3040 


1742-30 


1 


9 


6(3) 


4(26) 


7(15) 


13(4) 


81(32) 


— 


— 


— 


1 1 


1 


3 


65(140) 


1752-2806 


1749-28 


1 


.2 


3(4) 


-110(47) 


110(48) 


1.1(23) 


44(27) 


-4(6) 


-5(5) 


3 


0.7 


1 


.9 


630(280) 


1757-2421 


1754-24 


4 


.4 


-15(12) 


-520(760) 


520(760) 


-20(6) 


-500(400) 


— 


— 


— 


0.3 





.0 


11000(11000) 


1803 — 2137 


1800 — 21 


3 


9 


18(30) 


400(470) 


400(470) 


















7400(5700) 


1807-0847 


1804-08 


2 


.7 


7(9) 


1(19) 


7(10) 


-2.3(12) 


-4(5) 


-5(4) 


1(4) 


3 


0.2 





.0 


97(120) 


1818-1422 


1815-14 


7 


2 


-29(51) 


-87(210) 


92(202) 


-6(11) 


-78(65) 


— 


— 


— 


0.4 





.0 


3100(6800) 


1820-0427 


1818-04 


1 


9 


-14(2) 


26(7) 


29(6) 


-10(3) 


9(9) 


3(3) 


27(3) 


4 


3.1 





.1 


270(57) 


1824-1945 


1821-19 


4 


.7 


-12(14) 


-100(220) 


100(220) 


— 


— 














2300(3300) 


1825 — 0935 


1822 — 09 





.9 


11(23) 


39( 180) 


41(170) 






— 13(11) 


-9(5) 


3 


0.7 


o 


3 


170(710) 


1829-1751 


1826-17 


4 


.7 


22(13) 


-150(130) 


160(130) 


— 


— 














3500(2900) 


1832-0827 


1829-08 


4 


.9 


-8(9) 


16(20) 


18(18) 


-4(4) 


20(15) 


— 


— 


— 


0.3 





.1 


410(410) 


1833-0827 


1830-08 


4 


.7 


-38(26) 


-72(55) 


82(50) 


-37(4) 


0(15) 


— 


— 


— 


0.0 


1 


.0 


1800(1100) 


1835-1106 




2 


-8 


27(46) 


56(190) 


62(170) 


— 


— 


— 


— 


— 


— 






830(1600) 


1836— 1008 


1834— 10 


4 


5 


18(65) 


12(220) 


22( 140) 


















460(2900) 


1840 + 5640 


1839+56 


1 


.7 


-17(45) 


-29(43) 


34(44) 


-29(11) 


-33(12) 


-30(4) 


-21(2) 


1 


0.3 





.2 


270(350) 


1847-0402 


1844-04 


3 


.3 


1(15) 


5(38) 


5(38) 


-1(5) 


-9(19) 


— 


— 


— 


0.1 





.3 


83(590) 


1848-0123 


1845-01 


4 


.0 


-12(22) 


7(53) 


14(33) 


2(6) 


-41(19) 


— 


— 


— 


0.5 





.7 


270(590) 


1900-2600 


1857-26 


2 


.0 


-30(17) 


-35(67) 


47(52) 


-31(13) 


7(105) 


-19.9(3) 


-47(1) 


10 


0.6 





.2 


440(490) 


1913 — 0440 


1911 — 04 


2 


-8 


1(6) 


— 12(15) 


13(15) 


7(3) 


— 26(7) 


7(13) 


— 5(9) 


1 


0.7 


o 


,6 


160(200) 


1917+1353 


1915 + 13 


4 


.0 


0(12) 


-6(15) 


6(15) 


















110(280) 


1921 + 2153 


1919+21 


1 


1 


17(4) 


32(6) 


37(6) 


23(6) 


29(12) 


9(24) 


40(10) 


3 


0.6 





.2 


190(29) 


1932+1059 


1929+10 





4 


76(8) 


71(12) 


100(10) 


93.0(18) 


45(4) 


94.09(11) 


42.99(16) 


9 


2.2 


2 


.3 


180(17) 


1935+1616 


1933+16 


5 


6 


-6(2) 


-15(3) 


16(3) 


1.1(4) 


-16.3(7) 


-1(3) 


-13(3) 


3 


2.7 





.3 


440(87) 


1946 + 1805 


1944+ 1 7 


1 


.4 




2(10) 


2(10) 


7(1 ^1 


— 9(18) 


1 

L \ J ) 














1948 + 3540 


1946+35 


5 


.8 


-17(4) 


-7(4) 


19(4) 


-14(4) 


-3(4) 


-12.6(6) 


0.7(6) 


7 


1.0 


1 


.6 


510(76) 


1955 + 5059 


1953+50 


2 


.2 


-79(26) 


82(27) 


110(27) 


-30(3) 


59(30) 


-23(5) 


54(5) 


1 


1.8 





.8 


1200(280) 


2002+4050 


2000+40 


5 


9 


-9(21) 


-44(19) 


45(19) 


-17(9) 


-11(9) 


— 


— 


— 


0.3 


1 


.2 


1300(530) 


2013+3845 


2011 + 38 


8 


.4 


-97(22) 


-97(19) 


140(20) 


-37(9) 


-35(11) 


-32.1(17) 


-24.9(23) 


2 


2.8 


3 


.4 


5500(810) 


2018 + 2839 


2016+28 


1 


.0 


-13(12) 


1(3) 


13(12) 


-2.5(11) 


-4.6(15) 


2.6(2) 


6.2(4) 


8 


1.3 


1 


.7 


61(55) 


2022 + 2854 


2020+28 


2 


7 


-9(3) 


-15(5) 


17(4) 






-4.4(5) 


-23.6(3) 


8 


1.3 


1 


8 


220(54) 


2022+5154 


2021 + 51 


2 


.0 


-12(2) 


14(2) 


18(2) 


-5.7(10) 


-10.7(10) 


-5.23(17) 


11.5(3) 


8 


2.7 


1 


.2 


1700(17) 


2048-1616 


2045-16 





6 


54(41) 


220(150) 


230(150) 


141(14) 


-105(50) 


117(5) 


-5(5) 


3 


1.4 


1 


.5 


600(380) 


2108+4441 


2106+44 


5 


.0 


4(17) 


10(17) 


11(17) 


8(12) 


-10(11) 








0.1 





.7 


260(400) 


2113+4644 


2111 + 46 


4 


.5 


42(21) 


14(20) 


45(21) 


9(15) 


-1(16) 








0.9 





.4 


960(450) 


2157+4017 


2154+40 


3 


8 


35(19) 


11(19) 


36(19) 


0(17) 


7(15) 


17.8(8) 


2.8(10) 


2 


1.0 





.1 


650(330) 


2219+4754 


2217+47 


2 


2 


-28(8) 


-15(8) 


32(8) 


-14(3) 


-18(3) 


-12(8) 


-30(6) 


4 


1.3 





.3 


335(82) 


2257+5909 


2255 + 58 


4 


.5 


24(6) 


-7(5) 


25(6) 


28(7) 


-8(6) 








0.3 





.0 


550(120) 


2313+4253 


2310+42 


1 


.3 


33(6) 


11(6) 


35(6) 


21.2(15) 


6.3(15) 


-42(30) 


21(29) 


5 


1.5 





.7 


210(37) 


2321+6024 


2319+60 


3 





-7(30) 


4(30) 


8(30) 


-17(22) 


-7(19) 








0.2 







110(430) 


2326+6113 




4 




















0.4 






530(130) 


2354 + 6155 


2351 + 61 




4 


30(14) 


1(14) 


30(14) 


18(8) 


-1(6) 






1 


0.5 






490(220) 



References to previously published values are: (1) Harrison, Lyne & Anderson, 1993; (2) Brisken, 2001; (3) Fomalont et al., 1997; (4) Lyne, Anderson & Salter, 1982; 
(5) Fomalont ct al., 1992; (6) Han Kl Tian, 1999; (7) Brisken ct al. (2003); (8) Brisken ct al., 2002; (9) Chattcrjec et al., 2004; (10) Fomalont ct al., 1999. 



tempo standard deviation, versus the As parameter. Typ- 
ically, the range of As represented by the 74 pulsars in the 
sample is between —3.8 and 0, and the mean value of As is 
—2.2. The precise shape of the error factor curve depends 
on the whitening procedure, the power-law slope, the typical 



TOA uncertainties and systematic effects that may exist in 
the real data. 

To confirm the results of these simulations, we anal- 
ysed data for two pulsars with long data spans supplied to 
us from Jodrell Bank Observatory. Details of the data acqui- 



5 



(a) 







1- + 











10 

Time span (yr) 



15 



20 




-0.04 



(b) 



* 



10 

Time span (yr) 



15 



20 



Figure 2. Distribution of offsets in right ascension obtained from a timing fit to 100 simulated data sets dominated by (a) white noise 
and (b) red noise for time spans of 2, 3, 5, 10 and 20 years. The error bars are the mean formal error given by TEMPO for each of the 
data spans; see text. 



3 r 
2.5 

o 

O p . 

2 
c 
o 



1 - ; + 

0.5 I 1 1 1 1 1 1 

-3 -2.8 -2.6 -2.4 -2.2 -2 -1.8 -1.6 

Stability A 8 

Figure 3. The multiplicative factor required for the TEMPO un- 
certainties to equal the true uncertainties for different amounts 
of timing noise. Each point is obtained from 100 simulations of 
5 year data sets with TOA uncertainties of 0.5 ms. 

sition and analysis systems are given by Hobbs et al. (2004) . 
We first analyse a 20-year data set for PSR B0320+39 which 
contains little or no red timing noise. We select different data 
spans that are centred on MJD 49290 and fit for the pulsar 
position, setting the proper motion to the interferometric 
value. Right ascension offsets relative to the value from the 
full data span and their uncertainties as estimated by tempo 
are shown in Figure^. Results of a similar analysis using a 
35-year data set for PSR B1133+16, which shows a moderate 
amount of timing noise, are shown in Figure |3p. In line with 
common practice, higher-order frequency-derivative terms 
were used to whiten (i.e., absorb the red noise) from the 
final residuals when fitting for the position. Figure con- 
firms that tempo standard deviations are realistic error es- 
timates for a pulsar with little or no red timing noise for 
either short or long data spans. Figure 0p further illustrates 
that for noisy pulsars, data spans in excess of about 20 years 
are necessary for reliable position estimates and that tempo 
errors underestimate the true uncertainty by about a factor 
of two for data spans of the order of five years. 



Consistent with these results, for the Nanshan timing 
positions in Table 1 we have adopted an uncertainty of twice 
the formal tempo standard deviation for pulsars that show 
little timing noise (As < —1.8), whereas for pulsars with 
Ag > —1.8 we multiply the formal errors by three. 



4.2 Comparison of derived proper motions 

Values of e a and es given in Table 2 show that, with a few 
exceptions, the proper motions derived by comparison of 
Nanshan timing positions with previously catalogued posi- 
tions agree well with the best independently derived proper 
motions. This is further illustrated in Figure |S] which shows 
the distribution of (a) e a and (b) eg values. These plots 
exclude PSR B0736— 40; see below for a discussion of this 
pulsar. These distributions are close to Gaussian - 63% of 
the derived proper motions in right ascension lie within la, 
92% within 2a and 99% within 3a of the independent results 
- showing that error estimates are realistic. Similar results 
are obtained for declination where 63% are within lcr, 92% 
within 2a and 97% within 3a. 

However, there are a few conspicuous outlying points. 
For PSR B0736— 40, the scaled offsets e a and eg are very 
large, 11.9 and 5.5 respectively (Table 2). To derive the 
proper motion, the Nanshan timing position was compared 
with the timing position given by Downs & Reichley (1983) 
to give (i a = —48(2) and (i$ = 35(2) mas yr -1 . The 
e values are based on comparison of these proper mo- 
tions with the Brisken et al. (2003) interferometric proper 
motions given in Table 2: [i a — —14.0(12) and /is = 
13(2) mas yr -1 . Downs & Reichley (1983) independently 
determined a proper motion from their timing data: jj, a — 
—56(9) and fi$ — 46(8) mas yr -1 , and Siegman, Manch- 
ester & Durdin (1993) compared their Molonglo timing po- 
sition with the Downs & Reichley (1983) position to obtain 
fi a = -60(10) and fj, s = 40(10) mas yr -1 (Table 2). For 
this pulsar it is possible to solve for a proper motion using 
the Nanshan timing data alone, giving /j, a = —51(22) and 
fis — 6(24) mas yr -1 . Although not fully independent, all of 
these timing-based proper motions are in reasonable agree- 
ment with each other. An independent interferometric posi- 
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Figure 4. Right ascension offsets and their errors estimated using TEMPO for data sets of different lengths for (a) PSR B0320+39 and 
(b) PSR B1133+16. Data are from Jodrell Bank Observatory. 
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Figure 5. Distributions of offsets of derived proper motions from the best independently measured values relative to the estimated 
uncertainties (a) in right ascension and (b) in declination. PSR B0736— 40 has anomalously large offsets and is excluded from the plots. 



tion was obtained by Fomalont et al. (1997): fi a = —57(7) 
and us = —21(11) mas yr -1 . In right ascension this mea- 
surement agrees well with the timing results but it disagrees 
with the Brisken et al. (2003) value. In declination, it is 
discrepant with all previous measurements. This pulsar is 
near the southern horizon for northern interferometer sys- 
tems and corrections for ionospheric refraction are large and 
rather uncertain. We conclude that the timing results are 
probably more reliable, but this remains to be confirmed by 
future measurements. 

The proper motion estimates in right ascension and 
declination for PSR B2011+38 are also discrepant with pre- 
vious work. With the large estimated distance to this pulsar 
of over 8kpc our proper motion indicates an anomalously 
high transverse velocity > 4000 km s -1 . If we compare our 
position determination with that obtained by Hobbs et al. 
(2004) we obtain fx a = -28(17) and fj, s = -30(15) masyr" 1 
using the shorter 7-year baseline, but more precise posi- 
tion estimate. These values are consistent with previous 
estimates of the pulsar's proper motion and we therefore 



conclude that the uncertainty in the catalogue position was 
significantly underestimated. 



4.3 Pulsar velocities 

The two-dimensional or transverse velocity is given by Vr = 
4.74/itot-D km s _1 , where /Mot is the total proper motion 
in mas yr" 1 and D the distance in kpc. Figure |H] shows a 
histogram of the derived velocities (listed in Table 2) , omit- 
ting the anomalously high velocity for PSR B2011+38 men- 
tioned above and values which have an uncertainty of more 
than 2000 km s" 1 . This histogram has a similar form to ear- 
lier work. The mean and rms velocities of the sample are 
443 km s _1 and 224 km s _1 , similar to those given by Lyne 
& Lorimer (1994) (although these authors used the earlier 
Taylor & Cordes electron density model) and the new results 
by Hobbs et al. (2005) . The velocity distribution of our sam- 
ple favours a single-component Gaussian velocity model, as 
displayed by the dotted line in Figure||j] which is compatible 
with that of Lyne & Lorimer (1994) and Hobbs et al. (2005). 
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Figure 6. Distribution of pulsar transverse velocities for the pul- 
sars in our sample (values of the pulsar transverse velocity that 
have an uncertainty of less than 2000 km s ~ 1 are included) . 
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Figure 7. Projected directions of the proper motions for the 74 
pulsars in our sample relative to the Galactic Plane. The tail 
end of the vectors corresponds to the positions given in Table 1 
and the length of the vector represents the total proper motion. 
The effects of Galactic rotation and the unknown pulsar radial 
velocities have not been taken into account. 



PSR J1721— 3532 has a very high indicated velocity 
(7200 km s^ 1 ). While this is only a 3<r result, such a high 
velocity would be unique were it real. However, this result 
is more likely to result from a significant under-estimate 
of position uncertainties or an over-estimate of the pulsar 
distance. 

The motion of these pulsars in the Galactic plane 
is shown in Figure |7| where the length of the tracks 
indicate the pulsar velocities. This figure confirms that 
most pulsars are migrating from the Galactic plane 
JGunn fc Ostriker 1970t and only a few of them are mov- 
ing towards the plane jrlarrison, Lyne fc Anderson 1993| 
|Lyne, Anderson fc Salterl982*f ! However, for the major- 
ity of these pulsars, more precise proper motions al- 
ready exist in the literature. We also note that to make 
a full study of the pulsar space velocities, it is essen- 
tial to consider selection effect s (jHans en fc Phinney 1997| 
|Arzoumanian, Uhernoff fc Cordes 2002| . For instance, the 
Nanshan telescope only observed relatively bright pulsars 
in the the last few years and so high-velocity pulsars at 
large z- distances may have been missed (ICaraveo 19931 
|Frail, Goss fc Whiteoak 1994| . We therefore do not discuss 
the velocity distribution in detail and conclude this Section 
by discussing the new and improved proper motions from 
our sample. 



4.4 New and improved proper motions 

In Table |3] we list the nine pulsars for which we provide new 
proper motions (upper part of the table) and a further seven 
where our results are more precise than previous measure- 
ments. Besides the proper motions and distances repeated 
from Table 2, the 1-D velocities determined from Vm = 
4.74/xZ? where fi is the proper motion in right ascension or 
declination, their characteristic age r c = P/(2P) and surface 
dipole magnetic field strength B a = 3.2 x 10 19 (PP) 1/2 G 
are listed. The velocity uncertainties are determined solely 
from the uncertainties in the proper motion estimates and 



do not include any contribution from the distance estimates. 
For PSRs B0736-40 and B0844-35 these distances have de- 
creased by factors of 4.2 and 2.6 respectively with the new 
NE2001 distance model compared to the Taylor & Cordes 
(1993) model, significantly reducing the velocity estimates 
for these pulsars. Some of the derived velocities are very 
large, but in all cases, they are also very uncertain. Proper 
motion measurements for these relatively distant pulsars are 
difficult and more realistic velocity estimates await improved 
measurements. 



5 CONCLUSION 

We have presented new position and proper motion determi- 
nations for 74 pulsars based on five years of timing observa- 
tions at the Nanshan radio telescope of Urumqi Observatory. 
New or improved measurements are presented for 16 of these 
pulsars. It has long been known that the formal uncertain- 
ties given by the standard pulsar timing package tempo are 
often underestimated. Using realistic simulations of pulsar 
timing noise we show that the uncertainties in position es- 
timates given by tempo should typically be increased by a 
factor of 2 - 3. 

The Nanshan radio telescope is continuing to observe 
these and other pulsars. Over the next few years the uncer- 
tainties in the positions and hence the errors in their proper 
motions should be significantly reduced. The data set will 
also prove valuable in the study of glitches and timing noise. 
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Table 3. New and improved pulsar proper motions. 
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